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Abstract 

We compute the one-loop gauge couplings in six-dimensional non-Abelian gauge the- 
ories on the T 2 /Z2 orbifold with general GUT breaking boundary conditions. For 
concreteness, we apply the obtained general formulae to the gauge coupling running 
in a 6D SO (10) orbifold GUT where the GUT group is broken down to the stan- 
dard model gauge group up to an extra £7(1). We find that the one-loop corrections 
depend on the parity matrices encoding the orbifold boundary conditions as well as 
the volume and shape moduli of extra dimensions. When the £7(1) is broken by the 
VEV of bulk singlets, the accompanying extra color triplets also affect the unification 
of the gauge couplings. In this case, the B — L breaking scale compatible with the 
gauge coupling unification is sensitive to the change of the compactification scales. 
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1 Introduction 



Grand Unified Theories (GUTs) [1] have been reconsidered recently in the context of orb- 
ifold GUTs [2, 3] where orbifold boundary conditions in extra dimensions are utilized to 
break down a GUT gauge symmetry to the Standard Model(SM) gauge group and at the 
same time solve the doublet-triplet splitting problem. In orbifold GUTs, on top of the 
usual 4D logarithmic running to generate the difference between the SM gauge couplings 
at low energy [4], the Kaluza-Klein(KK) massive modes of the 4D gauge bosons give rise 
to additional threshold corrections. In orbifold GUTs, however, there is an ambiguity due 
to the existence of the non-universal gauge kinetic terms localized at the fixed points [5-7] 
where the local gauge symmetry is reduced compared to the bulk one. Nonetheless, by 
making a strong coupling assumption at the cutoff scale [8], the brane- localized gauge cou- 
plings may be ignored compared to the bulk gauge coupling, due to the volume suppression 
of extra dimensions. So, the orbifold GUTs can provide a minimal setup for considering 
the threshold corrections consistent with a successful gauge coupling unification. 

Over the past years, the 5D orbifold GUTs have been much studied as a simplest case, 
in particular, to compute the resulting effects of the KK massive modes to the gauge 
coupling unification for one flat [9] or warped [10] extra dimension. For larger GUT groups 
such as S'O(IO), however, the 5D case turns out not to be a minimal setup, because a usual 
Higgs mechanism is required for a further breaking to the SM gauge group or unwanted 
massless modes of extra components of gauge fields must get massive [11]. In contrast, the 
6D case has drawn more attention because it is more economic to obtain the SM gauge 
group directly from a large GUT group [12, 13] and there is more freedom to locate the SM 
fields for satisfying the experimental requirements such as the flavor structure [14] and the 
proton lifetime [15]. We may even regard the 6D orbifold GUTs as an effective theory of 
describing (heterotic) string compactifications as an intermediate GUT [16] with the hope 
to identify the remnants of string theory within the 6D orbifold field theory. 

In this paper, we consider the one-loop effective action for the gauge fields containing 
zero modes in six-dimensional M = 1 supersymmetric GUTs compactified on a T 2 /Z 2 
orbifold. This can be regarded as a generalization of the previous findings on orbifolds 
without gauge symmetry breaking [6] . In the presence of the orbifold boundary conditions 
that are commuting, we obtain the bulk and brane contributions due to bulk vector and 
hyper multiplets and identify the necessary counterterms to cancel the divergences appear- 
ing in dimensional regularization. A bulk vector multiplet leads to both brane and bulk 
corrections while a bulk hyper multiplet gives rise only to a bulk correction. The bulk 
divergences are cancelled by a higher derivative term with universal coefficient whereas the 
brane divergences are cancelled by brane-localized gauge kinetic terms the coefficients of 
which depend on the local gauge symmetry at the fixed points. In the case of the cutoff 
regularization [6] , there would be also power-like corrections in the cutoff scale to the gauge 
couplings, but they don't affect the gauge coupling unification at all. From the obtained 
effective action, we also derive the general expressions for the running of the effective gauge 
couplings for zero- mode gauge bosons. In the low energy limit, we consider the running of 
the gauge couplings, including the non-universal threshold corrections due to KK massive 
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mode£]. 

We apply the general formulae for the gauge coupling running in the six-dimensional 
50(10) orbifold GUT model proposed in Ref. [12]. This is the minimal setup to break 
50(10) down to the SM gauge group up to a U(l) factor only by orbifold boundary con- 
ditions without obtaining massless modes from the extra components of gauge bosons. In 
some realistic 50(10) orbifold GUT models, we discuss about the possibility of having a 
large volume of extra dimensions compatible with the success of the gauge coupling unifica- 
tion. We assume the breaking scale of the extra U(l) to be lower than the compactification 
scale in order to ignore the effect of the brane-localized U(l) breaking mass terms. 

For the case with isotropic compactification of extra dimensions, we show that the 
volume dependent term of the KK threshold correction can give a sizable contribution to 
the differential running of the gauge couplings for the large volume of extra dimensions. In 
this case, in order for the additional contribution due to extra color triplets to be cancelled 
by the volume dependent part, the breaking scale of the extra U(l) tends to be close to 
the compactification scale for the gauge coupling unification. On the other hand, in the 
case with anisotropic compactification, e.g. in the 5D limit where the bulk gauge group 
becomes the Pati-Salam SU(4) x SU{2) L x SU(2) R , we show that the shape dependent 
term of the KK threshold correction can be dominant, giving rise to the 5D power-like 
threshold corrections with non-universal coefficient in the compactification scales. These 
power-like corrections in the 5D limit are calculable, in contrast to the uncalculable power- 
like corrections in the cutoff scale in the genuine 5D case. Consequently, we show that 
the allowed contribution of extra color triplets or the breaking scale of the extra U{1) is 
sensitive to the shape modulus in a phenomenologically successful 50(10) model. 

The paper is organized as follows. First we give a brief review on the general boundary 
conditions for breaking the bulk gauge symmetry on T 2 /Z 2 . In Section 3, we present the 
one-loop effective action for gauge bosons in the general 6D orbifold GUTs and derive the 
running for the effective gauge couplings at low energy. Then, in Section 4, we consider 
the case with 50(10) bulk group and discuss the gauge coupling unification for some 
embeddings of the MSSM. Finally the conclusion is drawn. The details on the propagators 
on GUT orbifolds, the KK summations, the definition of special functions and some 50(10) 
group theory facts are given in the appendices. 

2 Boundary conditions on GUT orbifolds 

Before considering particular models, we give a brief sketch for the orbifold breaking of 
gauge symmetry in a six- dimensional non-Abelian gauge theory with a simple gauge group. 
Two extra dimensions are compactified on the orbifold T 2 /Z 2 . For the extra coordinates 
z = x 5 + ix 6 , there are double periodicities z ~ z + 2n(R 5 n 5 + iRqUq) with radii R 5 , R 6 and 
integer numbers n 5 ,n 6 . Further, when the bulk positions are identified by a Z 2 reflection 
symmetry as z — > —z, there are four fixed points on the orbifold: zq = 0, z\ — ttR*,, 
z 2 = inRs and z 3 = ttR 5 + ittRq. 

1 For some early works on string theory computation of the one-loop gauge couplings, see Ref. [17, 18] 
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In order to break the bulk gauge symmetry down to the SM gauge group, let us in- 
troduce nontrivial boundary conditions for bulk gauge fields A M with M — 0, 1,2,3 = \i 
and M = 5, 6 = m. The boundary conditions are specified by unitary parity matrices 
Pi(i — 0, 1, 2, 3) at the fixed points, 



ij ■ 



ViA^z)Vr = PtA^-z + Zi )Pr l = A^z + z, 
TVU^V 1 = -PiAmi-z + z^Pr 1 = A m (z + Zi ) (1) 

where Pf — l(i — 0,1,2,3). The above boundary conditions can be rewritten simply in 
terms of component fields with Am — A a M T a in the group space as 

Aft-z + Zi) = {Q l ) a b A b ^(z + z^, (2) 
A a m (-z + Zi ) = -(Q t y h A b m (z + Zl ) (3) 

where 

(Qi) a b = tT(T a P i T b P i ). (4) 
Here the defined matrices {Qi) a &(i = 0, 1, 2, 3) fulfill 

(Qi) a a'(Qi) b b'Vab = Va'b', fabc(Qi) a a' (Qi) V (Qi)° c' = fa'V c' (5) 

where r) a t, is the Killing metric defined by tr (T TJ,) = r]ab on the group space and it is used 
to raise and low adjoint indices, and f a b c are the group structure constants given in the 
group algebra [T a ,T 6 ] = if a b c T c . Note that Q\ = 1 from the Z 2 symmetry and hence Qi 
are real symmetric matrices. Eq. (J3J) and the second property in eq. §5§ can be rewritten, 
respectively, as 

P t T a P t = {Qi) a b T\ (6) 

Q l T£Q l = (Q.nr* (7) 

with (T b ) ac = if abc . 

We also discuss on the Wilson lines on a torus in comparison to the local boundary 
conditions as given above. The boundary conditions along noncontractible loops on a torus 
are defined by the unitary matrices U\, 11% as 

U x A M {z + 2nR 5 )U^ = A M {z), (8) 
U 2 A m (z + i27tR 6 )U^ = A M (z). (9) 

Since x 5 + ttR 5 — > —x 5 + nR 5 is equivalent to x 5 + ttR 5 —> —x 5 — tcR 5 —> — x 5 + ttR 5 and 
similarly for the other coordinate, we obtain the following relations, 

U 1 = P 1 P , U 2 = P 2 P . (10) 

Then, we can see that the consistency conditions for the Wilson lines on orbifolds, U\PqUi = 
Po and U 2 PqU 2 = Po, are satisfied. Moreover, since U 2 UiPq = P3 and [Ui,U 2 ] = 0, the 
parity matrix P3 can be written as 

P 3 = P 2 P P! = PlP P 2 . (11) 
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Therefore, the Wilson lines one can consider are not independent of local boundary condi- 
tions, and one of the parity actions is not independent. 

For simplicity, let us focus on the case with commuting parity matrices, i.e. [Pi, P 3 ] = 
or [Qi, Qj] = 0. For these parity actions, the rank of the gauge group is not reduced. In this 
case, it is convenient to choose the Cartan-Weyl basis such that the orbifold actions become 
diagonal. In this basis, the generators are organized into Cartan subalgebra generators 
Hi, I = 1, • • • , rank(G), and the remaining generators, E a , a = 1, ■ ■ ■ , (dim(G) — rank(G)), 
with 

[Hj, E a ] = atjE a , (12) 

where ai is the rank(G)— dimensional root vector associated with E a . Then, it is always 
possible to write the parity matrices as 

p. = e -am*.fl ( 13 ) 

which defines the rank(G)— dimensional twist vector Vi for each fixed point. Thus, the 
relations ([6]) become 

PHjP = Hj, (14) 

PE a Pi = e~ 2 ^E a . (15) 

In this basis, the matrices Qi are also diagonal such that (Qi) 1 j = 6j and (Qi) a p = 
e -2ma-Vi fia anc j Q^jjgj. entries are zero. Here we have that a ■ Vi — or | mod Z for Z<i 
actions at the fixed points because Qf — 1. Then, a bulk field takes a combination of 
parity eigenvalues (po,Pi,P2) with pi = +1 or —1 under three independent Z 2 actions, so 
it is composed of a subset of basis functions on a torus with radii 2i? 5 and 2i? 6 . 



3 The effective action on GUT orbifolds 

In this section, we present the general formulae for the one-loop effective action in a 6D 
M = 1 supersymmetric GUT where the bulk gauge symmetry is broken by local boundary 
conditions at the fixed points on the T 2 /Z 2 orbifold as described in the previous section. 
As a result, we also discuss about the running of the gauge couplings for zero- mode gauge 
bosons at low energy. 

3.1 The one-loop effective action on the T 2 jZ<i orbifold 

We consider a 6D M = 1 supersymmetric non-Abelian gauge theory compactified on the 
orbifold T 2 /Z 2 . In terms of component fields, a vector multiplet is composed of gauge 
bosons A M and (right-handed) symplectic Majorana gauginos A while a hyper multiplet is 
composed of two complex hyperscalars <p± without opposite charges and a (left-handed) hy- 
perino ip. Since all charged hyperinos have the equal 6D chiralities due to supersymmetry, 
one is not allowed to write the 6D mass terms for hyper multiplets. 
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In the process of taking the usual gauge fixing for a non-Abelian gauge theory [6], we 
also introduce ghost fields c a . Then, the orbifold boundary conditions for bulk component 
fields that we are considering are as the following, 

Afa,-z + Zi) = (Q l ) a b A b l (x,z + z i ), A a m (x,-z + z i ) = -{Q i ) a b A b m (x,z + z i ), 

c a (x,-z + Zi ) = (Q t ) a b c b (x,z + z t ), X a (x,-z + z t )=i 1 5 (Q i ) a b X b (x,z + z t ), 

ip(x,-z + Zi) = ry 5 rjiPii(j(x,z + Zi), 

(f>+(x,-z + Zi) = rjiPi(j) + (x,z + z^, cj>-(x,-z + = -r)i<f>-(x,z + Zi)Pi (16) 

with i — 0,1, 2, 3. Here the forms of the parity matrices depend on the representation of a 
hyper multiplet under the bulk gauge group. Each hyper multiplet can take its own value 
of r/i as either +1 or —1. 

Taking into account the group structure of propagators in loops as discussed in the 
Appendix A and following the similar procedure as in the case with no orbifold breaking of 
the gauge symmetry in Ref. [6], we obtain the one- loop effective action for the background 
gauge bosons up to quadratic orders as 

r<2) [4J = E / j0ji A X- k - A "<- k < H~ <* 2 " + *"*") 



k,k' 



(17) 



x<^ -n^ + 4(A;V-^) n ++ - 2k ■ k'g^(U%_ + U 



tG 



where U% = + + n*,, ng, = U h ^ + + U h ^_ + Tl% with 



{^u,±) a b 



= E 



p,p 



d A p 



Tr 



- (2p + k)^(2p + k) v G g ,±{p + kj + k' , p + k) 



(18) 



Ku,±)\ = E/ J^ Tt {-(2p + kU2 P + k)„G h>± (p + kJ + k>,p + k) 



p,p 



+ 2i 9^ 6 pi,p+k-k>} T bGhAP, Pi f) T ° 



p,p 



d A p 



Tr 



D^(p,pJh^T b D4p + kJ + k',p + k)-i u T a \, (21) 



(19) 



/?4 
— ^Tr [b x {p, p, i?h,T b D x (p + k,ff + k', p + k) lu T a ] , (20) 



and 



E 



p,p' 



d A p 



G g ,±(P + k,P + k,ff + k')T b G g , ± (p,pJ)T a 
G g , T (p + k,p + kj + k')T b G g , ± (p,pJ)T a 



(22) 



= (IT 



±,±) b- 



(23) 



Here the propagators appearing in the loops are given in the Appendix A. Since we consider 
the commuting parity matrices, the orbifold actions with respect to the fixed points other 
that the origin are factorized out of the propagators which would be given for the case 
with one Z 2 orbifold action only. Then, after identifying various equivalent terms, we get 
the effective action in a simpler form as a decomposition into bulk and brane parts, 



bulk + r 

branc 



(24) 



with 



bulk 



k,k' 



X 



-^5 a b -i(U G + U H )\(k,k) 
d A k 



branc 



2E J 7^4(-^-^)^(^^)(^-^)[-4i(n G )%] 

k,k' 



(25) 
(26) 



where 



(u G ) a b (k, k) = ^ y [ 

{^) d {p 2 -f)[{p+kf-{p+kf} 

x itr Adj [{ 1 + cos(2p 5 vr J R 5 )QoQi } { 1 + cos(2p 6 nR (i )Q Q 2 ^T a T b j , 

(iih)V*,£) = V£ / 70 1 



(27) 



x 



(27T) d (p2 _ p2) [(p + fc )2 _ (^ + £)2] 

itr R [{l + Wl cos(2p 5 7r J R 5 )P P 1 }{l + W2 cos(2p 6 7r J R 6 )P P 2 }T a T f) ], (28) 



(n G )%(M',£) 



/i 



4-<2 



E 



(^) d ( P 2-pP)[(p + ky-(p + ky} 



X 



itr Adj [{l + cos(2p 5 7r J R 5 )g Q 1 }{l + cos(2p 6 7iR 6 )Q Q 2 ^Q T a T b j . (29) 



Here /x is the renormalization scale in dimensional regularization with d = 4 — e, and 
P = (P5,P6) = ^) with 715, n6 being integer, and similarly for k and k! . In simplifying 



the expressions in the above, for the generators satisfying P{F a = ±T a Pj, we made use of 
cos(2(j» 5 + k 5 )irR 5 ) = ± cos(2p 5 7rP 5 ) and cos(2(p 6 + fc 6 )7rP 6 ) = ± cos(2p 6 7rP 6 ). Further, we 
notice that trAdj is the trace over indices of the adjoint representation and tr^ is the trace 
over indices of the R representation. 

From eq. (!24"1) . we can see that a vector multiplet gives rise to both bulk and brane- 
localized corrections while a hyper multiplet only leads to a bulk correction. It has been 
shown that the absence of the brane-localized corrections due to a hyper multiplet is 
restricted to the case with even ordered orbifolds [5]. 

In order to simplify the expression for the bulk contribution (1251) . we define the quantity 

n (p 5 ,p 6) _ J dd p 1 



i 



(27r) d (p2_p2)[( p+ fc)2_(£ + £)2] 

(27r/i) e / dxJo[x{l-x){k 2 + k 2 ),xk 5 R 5 + p 5 ,xk 6 R 6 + pQ\ (30) 
(4tt) z V J 

where p = ( n5 ^f s , rab ^" 6 P6 ) with p 5 , p 6 = or \ and n 5 , n 6 being integer, V = (2n) 2 R 5 R e , and 

e -TTt[c+a 1 (n 1 +c 1 ) 2 +a2(n 2 +C2) 2 } /g-g 



/CO 



-e/2 



with a>i = l/R 2 +4 (i = 1,2). Thus, we can rewrite the bulk contribution due to a vector 
multiplet as 



G b 



hr Adi [T a T b ] (n<°'°> + n^) + n&°> + n (hh)) 
+\ti Adi [QoQiT a T b ] (n(°'°) + n^D - n(!>°) - n<H>) 
+\tr Adi [Q Q 2 T a T b ] (n^°) - n^) + n(i'°) - n<™>) 
+|tr Adj [(5iQ 2 T T 6 ] (n(°'°) - n(°-5) - n&°> + n&*>). (32) 

Similarly, we can write the bulk contribution of a hyper multiplet as 

(n H ) a b = ~tr R [T a T b ] (n(°'°) + n(°-5) + n&°> + n.(U)) 

-\wn tT R [P PiT a T b ] (n(°'°) + nf°'3) - n^°) _ n^>) 
-J w2tr R [p p 2 T o r 6 ] (n(°'°) - n< ^ + n(i°) - n&*>) 
-^x^trR [Pi p 2 T a T b ] (n (0 ' 0) - n (0 '5) - n(s.o) + n^>§)). (33) 



Also defining 



n(p 5 ,P6) = HiLly* [ _fl 6 -2p,k-k' (M) 

2 (27T) d ( p 2 _p2)[(p + A;) 2 - (p + £) 2 ] 



S 



where p = ( " B r 5 P5 > " 6 ^" 6 P6 ) with ps,p6 = or | and n 5 ,n 6 being integer, we can rewrite the 
brane contribution as 

(n G )% = ^tr Adj [g T a T 6 ](n( ' ) + n( ^) + n(i ) + n(i^) 
+^ti Adi [Q 1 T a T b ] (n(°'°> + n (0 ^ } - n^°) - nd.D) 
+^tr Adj [Q 2 T a T 6 ] (n(°'°) - n(°'i) + n^°) - n<™>) 
+^tr A dj[Q3T a T 6 ] (n(°>°) - n(°-3) - n(3-°) + n&s)). (35) 

Thus, we find that the brane-localized contribution at each fixed point corresponds to the 
brane projection of the bulk quantity by the local parity matrix. 

3.2 Counterterms for loop divergences 

After the KK summation given in the Appendix B, we can separate the divergent term as 



TT(P5.P6) 



7T 



(An) 2 V 6 

Thus, the bulk contribution becomes 



-R 5 R 6 (k 2 + k ) 



-2 



+ O(e ). 



(n G + n H y b = htr Adi [T a T b ] - trR [T a T b ])—^—^R 5 R 6 (k 2 + k 2 

4 {47T y V O 



-2 



(36) 



+ O(e ). (37) 



Therefore, we find that the loop corrections generate a divergent higher derivative term 
the coefficient of which is proportional to the universal M = 2 beta function. At the 
momentum scale higher than the compactification scales, the higher derivative operator 
becomes important so that the gauge couplings run power-like in momentum scale rather 
than logarithmically [6] . However, it does not affect the unification of the gauge couplings, 
even if it is important in determining the value of the unified gauge coupling and the 
unification scale. 

For a given set of ingoing and outgoing momenta of gauge bosons satisfying p = k " ' 
we compute IT P5 ' P6 ) as 



2 ' 



d^{i +M4 ^ e "" ] -J! dxln 



x(l-x)(k 2 +k 2 ) + [—+(x — )k 



k' 



1 



(38) 



Thus, the brane contribution becomes 



EtrAdj^T^] 5 2 {z - Zl ) + O(e ). 



i=0 



(39) 



Therefore, the appearing divergent term at each fixed point respects the corresponding 
gauge symmetry which depends on the local orbifold action. 
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Consequently, in order to subtract the e poles in (Uq) 0, b, (J^h) 0, b and (JIg) 11 b obtained 
in eqs. (1371) and (|39|) . we require the following new counterterms which are not present in 
the original action, 



Cct. = J d 2 zd 2 e ^ti[wn 6 w] + (^-w a w a 8 2 ( z - Zi ) 



+ h.c. 



(40) 



Here h 2 is a dimensionless bulk coupling while g 2 a {i = 0,1,2,3) are dimensionless gauge 
couplings correspoding to the local gauge groups at the fixed points. Note that the brane- 
localized gauge coupling can be non-universal so it could affect the predictive power of the 
orbifold GUTs. 



3.3 Limiting cases 

For a later use in the running of the zero-mode gauge coupling, let us take the asymptotic 
limits of the loop corrections. First, in the low momentum limit k 2 <C 1/R 2 6 , with k = 



k' = 0, by using eq. (1B.4|) . we get the approximate forms for eq. fl30l . 



n(°>°) 



n(°>!) 



n(i-l) 



(Att) 2 V { 6 



7T 



RsRsk 2 



-In 



ne^ E fi 2 Rl\r](iu) 



In 



47rV 2 |r ? (m)| 4 i?^ 2 
2 



(Att) 2 V { 6 



7T 



R 5 R 6 k 2 



ln(7re 7B /i 2 ^) _ nu + 2 ^(-lfln |1 - e 

n>l 



(41) 



1-Kun 1 2 



-In 



#i(l/2|m) 



rj(iu 

7i 



(An) 2 V { 6 
-In 



R 5 Rek 2 



-2 



ln(4vre 7B /i 2 i? 2 ) -2 J^ln 



n>l 



1 + e 



1 - e" 



2n 



^(-m/2|m) c _ W4 
r;(iw) 

2 



7T 



(47r) 2 V i 6 
-In 



R 5 Rek 2 



ln(47re 7B /i 2 J R 2 ) - 2 J^(-l) n In 



n>l 



1 + e" 



1-e" 



A(l/2-m/2|m) c _ W4 



(42) 



(43) 



(44) 



Further, using; the fact that n(°>°) + II^) + nM) + n^-l) is the same as the KK sum on 
a torus with each radius double sized and with no Wilson lines, i.e. from the approximate 
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form of the sum, 



P5,P6=0,i 



(Att) 2 V { 6 



7T 



-(4R B R 6 )k 2 



e 



In 



167r^e-^|r7(i7x)I 4 ^A; 2 



(45) 



we note the useful identity for the theta functions, 



#i(l/2|m) tfi(-m/ 2 M #i(l/2 - iu/2\iu) = A\r](iu)\ b e 



I 6 „7YU 



(46) 



In the high momentum limit k 2 3> 1 / i?§ 6 , with k = k' = 0, eq. (1301) becomes, indepen- 
dently of the orbifold actions, 



JT(P5,P6) 



1 l-R 5 R 6 k 2 



(4n) 2 V \6 



e fc 2 



(47) 



Therefore, from eqs. fl32l) and fl33|) . even the finite part of the bulk correction becomes 
universal at high energy. 



3.4 Running of the 4D effective gauge coupling 



In this section, we consider the running of the effective gauge coupling which is defined as 
the coefficient of the kinetic term of a zero-mode gauge boson. It also includes the bulk 
higher derivative term and the brane kinetic terms. 

From the one-loop effective action (1241) . the zero- mode gauge coupling reads 



1 



9 'tree, ab 



k 2 V 

1 Hree 



6 ab + iV(U G (k, 0) + U H (k, 0)) a b + 4i(U G (k, 0, 0)) a b (48) 



with 



9tree,ab 



v_ 

~n2 



i=0 y *> a 



Jab- 



(49) 



M 2 



When taking the minimal subtraction scheme for divergences ( 1371) and ( 1391) at k 2 
where M* is the 6D fundamental scale, we define the renormalized bulk and brane gauge 
couplings for £i/z 2 = M 2 with £i = 47re 2_7B at that scale. Then, below the compactification 
scales (k 2 < l/i?| 6 ), using eqs. (J38J) , (pEI]) -([4l | with flU, we have eq. (gHJ) as 



9ett,ab(k' 2 



1 A 1 R 1 M * 1 >^ R« T 1 

W a + l6^ Bab M I 2 " " 16^ 2^. ~ 4^ a& 



(50) 



i,j=± 
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where g rA are the renormalized gauge couplings and the beta functions are 



B, 



ab 



- tr Adj (T a T b ) + ^tr R (T a T 6 ) - 2 ^ tr Adj (Q z T a T b 

R i=0 

, 3 



(51) 



and 



with 



B + ab + 


1 

4 


-tr Adj (T a T 6 ) + ^tr R (T a T 6 ) , 

R 


(52) 


B- ab + 


1 

4 


- tr Adj (Q QiT a T fe ) + ^r7 VtrR(P PiT a T 6 ) 

R 


(53) 




1 

4 


- tr Adi (Q Q 2 T a T b ) + J2voV?tr R (P P 2 T a T b ) 

R 


(54) 


B ab~ 


1 

4 


- tr Adj (QoQ 3 T a T 6 ) + ^2 V ^tr R (P P 3 T a T b ) 

R 


(55) 



L-+ 



In 
In 
In 
In 



4e- 2 \ri(iu)\ 4 uVM* 

2 1 4 



L 4 



tfi(-|m) 



7?i(-^m|m)e- ,ra/414 



Mi/Mr 



^(- - -m|m)e-'™/ 4 



uVMt 



(56) 
(57) 
(58) 
(59) 



Further, K ab corresponds to the power-like dependence on the momentum scale and it is 
suppressed by the compactification volume at low energy as in the case without orbifold 
breaking of the gauge symmetry [6]. B ab are the M = 1 beta function coefficients of the 
logarithmic running due to the massless modes. They are composed of both bulk and 
brane corrections. On the other hand, B a \ are the M = 2 beta function coefficients for 
the KK massive modes of the bulk fields. The logarithms Lij have the common volume 
{V) dependence, but also they are functions of the shape modulus (u), being of different 
form depending on the parities. Since the KK massive mode correction contains a non- 
universal part due to the gauge symmetry breaking, they can affect the unification of gauge 
couplings. 
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4 Gauge coupling unification in a 6D 5*0(10) orbifold 
GUT 



We consider a 6D J\f = 1 supersymmetric SO (10) orbifold GUT and compute the gauge 
coupling running by using the general formulae found in the previous section. 

4.1 Orbifold breaking of 50(10) 

In order to break the bulk SO (10) gauge group down to the SM one, we introduce the 
parity matrices in eq. (pQ) or (fl6l) for a fundamental representation [12] as 

Po = Wo, (60) 
P 1 = diag(-l,-l,-l,l,l)x<7°, (61) 
P 2 = diag(l,l,l,l,l)xa 2 , (62) 

and P3 = -P1-P2 from the consistency condition (jlip . Then, the parity operations Pi,P 2 
break SO (10) down to maximal subgroups, the Pati-Salam group S77(4) x SU (2) l X SU (2)r 
and the Georgi-Glashow group SU(5) x U(l)x, respectively. The parity operation P 3 also 
breaks SO (10) down to the flipped SU(5) but it is not an independent breaking. Thus, 
the intersection of two maximal subgroups leads to SU(3)c x SU(2)l x U(1)y x U(1)x as 
the remaining gauge group. This can be seen from the gauge bosons with positive parities: 
45 is decomposed into (15, 1, 1)+ + (6, 2, 2)_ + (1, 3, 1)+ + (1, 1, 3)+ under Pi (where ± 
indicate the parities) and 24 0j + + 10_4 5 _ + 10 4 _ + 1 ,+ under P 2 . Then, finally, the extra 
U(l)x or U(1)b-l has to be broken further by the VEV of bulk or brane Higgs fields. 

For applying the parity action to other representations, from eq. (fT3"j) . we can rewrite 
the parity matrices in terms of Cartan-Weyl generator^ as a special case of eq. (fT3"j) . 

P 1 = e -a^i(-«V+Tx) j ^ * (63) 
P 2 = e~ 2mX2T *, x 2 = l (64) 

o 

where Ty, Tx are the U (l)y and U(l)x generator^, respectively. We consider a set of hyper 
multiplets, N w 10's and N\e 16's satisfying iVi = 2 + N w for no irreducible anomalies 
[21,22]. Both N w and iV 16 have to be even for the absence of localized anomalies unless 
there are split multiplets at the fixed points [22]. A 10 = (H, G, H c , G c ) is decomposed into 
(6, 1, 1)_ + (1, 2, 2) + under P 1 and 5_ 2 ,- + 5 2 ,+ under P 2 . Then, we get a massless Higgs 
doublet from H c of 10. On the other hand, a 16 = (Q, L, U, E, D c , N c ) is decomposed 
into (4, 2, 1) + + (4, 1, 2)_ under P 1 and lO^- + 5_ 3j+ + l 5j+ under P 2 . Then, we also get 
a massless lepton doublet from L of 16. The M = 2 partner of each hyper multiplet has 
the parity matrices in the group space multiplied by the negative overall parity, due to the 
discrete choice of the Scherk-Schwarz twist in SU(2) R space. We note that in eq. (FlBl) . 
r]o = 1 and r/3 = for each hyper multiplet. 

2 Onc has to be careful with multiplying a U(l) phase for the correct parity matrices satisfying Pf = 1. 
3 See the appendix D for details. 
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4.2 Gauge coupling running at low energy 

In order to break the extra U(l)x gauge symmetry by a usual Higgs mechanism, one can 
introduce 16 Higgs multiplets in the bulk [13,14]. In that case, after the orbifolding, on 
top of SM singlets, one ends up with extra color triplets as zero modes. Since the extra 
color triplets can get masses of order the B — L breaking scale Mb-l at the fixed points, 
they already start contributing to the running of gauge couplings at that scale. Thus, from 
the general result ( l50i) . we consider the logarithmic running due to zero modes by taking 
two steps across the B — L breaking scale. 

The brane-localized B — L breaking masses for the color triplets can modify their KK 
massive modes so that there exists an additional contribution to the the gauge coupling 
running. However, when the B — L breaking scale is below the compactification scale, the 
new contribution becomes suppressed as M B _ L /M^ with M c = 1/y/V. Thus, henceforth 
we assume this case to ignore the effect of the brane-localized B — L breaking masses. Then, 
much below the compactification scale, the running of the 4D effective gauge coupling of 
the SM gauge group is governed by 



9cB,ab(k 2 



9a 



'ab 



1 , M B _ L 
7^> — 



+ 



1 s . Ml 



167T 



M 2 B _ L 



16tt 2 ^ 

h3=± 



B1U 3 + 



57T 



:AJ, 



a<Jab 



(65) 



where g u is the universal renormalized gauge couplings and A a are corrections due to renor- 
malized gauge couplings localized at the Pati-Salam and flipped SU (5) fixed points. B' ab = 
b a S a b are the beta functions of the gauge couplings in the MSSM with b a = (33/5, 1, —3) 
given below the B — L breaking scale. Moreover, above the B — L breaking scale, the beta 
functions for the gauge couplings are given by B a b = B ab — C a b + B a b where 



B, 



1 r 



ab 



-3tr Ad j(T a T 6 ) + ^tr R (T a T b 



R 



+7 Z) [ " 3tr Adj (Q*T a T 6 ) +^7 7 ftr R (P i T a T 6 )" 



(66) 



i=i 



R 



with Q 3 = Q1Q2, P3 = P1P2 and 77^ = r/fr?^ , and C a b = c a 5 a b is the contribution coming 
from vector-like massless modes which get tree-level brane masses of order the GUT scale 
and B a \j = b a 5 a b comes from the brane-localized fields. Further, the beta functions of the 



4 Although there are also power-like threshold corrections in the cutoff regularization [6,9], they don't 
contribute to the differential running of gauge couplings. Nevertheless, the power-like contributions may 
have the net effect of placing an upper limit on the possible volume of the extra dimensions [23] . 
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KK massive mode corrections are 





1 

4 


-tr Adj (T a T 5 )+^tr R (T a T 6 ) , 

R 


(67) 


B- b + 


1 

4 


- tr Adj (g 1 T a T 6 ) + ^7 7 ftr R (P 1 T a T 6 ) 
i? 


(68) 


B+ ah - 


1 

4 


- tr Adj (Q 2 T a T b ) + ^^tr R (P 2 T a T b ) 
i? 


(69) 


B ab 


1 

4 


- tr Adj (Q 3 T a T b ) + ^^tr R (P 3 T a T 6 ) 


(70) 



R 



4.2.1 Computation of traces 

Now we compute the necessary traces to get the running equations. To this, we define the 
following invariant quantity including all SO (10) gauge fields, 

D — D T? a T?b 

= ^ab^ SO{W) P SO{W) 

= B V + B M (71) 

where 

3 3 
B v = -- [tr Adj F| o{10) + ^tr Adj (QjF| o{10) )j , (72) 

i=i 

1 3 

Bm = Ei[ tr ^o(io) + E^ tr «( P ^o(io))]- (73) 

ft i=l 

Moreover, similarly we define 

^' = ^50(10)^0(10)- (74) 

By using the traces for maximal subgroups of 5*0(10) in the Appendix D, we obtain 
the following result for the vector multiplet, 

3 

Bv = -^tr Ad jF 5O(10) + Qtr 2 F su(2)L + 6tr 2 F S!7(2)j? 

-Str 5 Fg U(5) + 6F^ (1)x - Str 5 Fg U(5y + 6F u{1)xi . (75) 

For the hyper multiplets, we also have 

B M = B w + B 16 (76) 

15 



with 



Bio 



- T^iotrio-Fjo(io) + 9 X [ ~~ tr 4-^li/(4) + tr 2 -F| C /( 2 ) jL + tr 2 F| t/(2)jl 

10 

= ^^letrie-^Jocio) + X ^ [ tr2 ^5C/(2) L ~ tr 2 F su(2) R 

16 

+ lE^ 16 [- 2tr A(5) + ^(Dx" 
4 16 L 



Further, we have 

B ++ 
B- + 



^[-trAdj-^lo(io) + ^iotrio-Fjo(io) + Ni 6 tri 6 F| o(10) 
2tr 2 Fg U(2)L + 2tr 2 F| [/(2)jj 

+ 2 ^ [ ~ iliF SU(i) + tr 2-^lc/(2) L + tT 2^SU(2) R 
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+ XI ^ [ tr 2^ 7 'li7(2) i - tr 2-^lc/(2)j 



5" 



B~ 



16 



_tr 5-^!sC/(5) 



2 + 2 ^(D* 



? 2 

51/(5) 



172 

+ 2 uWx 



- J> 2 16 [ - 2tr 5 F, 

16 

-tr5-^5C/(5)' + 2 F u(i) x , 



16 



(77) 



(78) 



(79) 



SO) 



(82^ 



Therefore, by reading off the gauge kinetic terms for the SM gauge group in B and B 1 *, 
we can find the general expression for the beta function coming from the massless modes 
in the bulk as B ab = b a 5 a b with 

ba = b v a + b™ + (83) 



where 



lV 



(0,-6,-9), 

^ 10 (1, 1,1) + 1X^(1 1,-1), 



(84) 
(85) 
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bl 6 = ^(27V 16 -^^)(l,l,l) + I^^(-^2,0) 

16 16 



hE«([.-i.-D 

16 



(86) 
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and the beta function for KK massive modes as B^ b = b l j5 a b with 

b + a + = ~(-8 + iV 10 + 2iV 16 )(l,l,l), (87) 
K + = ^4 ) 0) + i^,f(il,-l) + i^,f(-|2,0) ) (88) 

10 16 



K~ = ^(2 + 5> 2 16 )(-l,-l,-l), (89) 

16 

K- = ^(f,-2,-2) + i^« 6 (I -1,-1). (90) 

16 

Here, in order to get the beta function for U(1)y, we made use of the relations between 
£7(1) gauge bosons (ID. 231) in the Appendix D. The corrections due to hyper multiplets in 
eq. (1771) and flTHj) also contain mixing terma^l between the U(l)y and U(\)x gauge bosons. 
After transforming to the canonical gauge kinetic terms, this mixing leads to an overall 
shift in the U(l)x charges as well as the coupled renormalization group equations for two 
U{1) gauge couplings and the U{l)x charge shift [19]. However, when the extra gauge 
boson gets a heavy mass for giving the see-saw scale for neutrino masses, the mixing effect 
is not relevant for the low energy physics while the running of the gauge coupling for a 
light U(l) gauge boson is not affected by the presence of the mixing term [19]. 

Consequently, compared to eq. (1831) . we obtain the relation between beta functions as 



b a = (0, -4, -6) + b + a + + b- + + 6+~ + b--. (91) 

The first term is only due to the difference between the beta functions of M = 1 vector 
multiplets and M = 2 vector multiplets for the SM gauge group. Apart from that, the sum 
of the M = 2 beta functions for the volume dependent part of the KK massive contributions, 
i.e. ^ijWj, contains only the KK massive modes for the bulk fields containing the zero 
modes. Therefore, from the beta functions (1851) . (1861) . (IHBl and ( 190|) . one can find that the 
terms proportional to r/f or r?f?7<f , i.e. the orbifold actions associated with Pati-Salam and 
flipped SU(5) gauge groups generate the non-universal corrections to the gauge couplings. 
From the obtained beta functions ( 183]) and (l871)-( l90j) . eq. ( |65|) becomes 

4tt 4t 1 ; , Ml 1, Mir 1 A, 



+ ^~b a In • 



Here b a = b a — c a + b a is the M = 1 beta function above the B — L breaking scale and b l J 
are the beta functions for the KK massive modes. 



5 The gauge kinetic terms localized at Pati-Salam and flipped SU(5) fixed points can also lead to a 
mixing. 
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4.2.2 The differential running of the gauge couplings 

For a number of hyper multiplets with arbitrary parities, we assume that both vector-like 
particles (getting brane masses of order the GUT scale) and brane-localized particles fill 
GUT multiplets, i.e. c a and b a are universal. Then, we get the general formula for the 
differential running of gauge couplings as 



12 1 5 1 



1 



where 



1 <lln M 



9 9 



M B _ L 2 + 2 



15 



A 



1 = ?-n£»:°-S5>; 6 ^£«> 

10 16 16 



9 9 



15 



1 - _ Vr7 10 - — V 

7 14 11 14 



10 



14 



16 



16 



16 



16 



(93) 

(94) 
(95) 
(96) 



Thus, we find a general relation between coefficients as 

b = ^ + b- + + b—. 



(97) 



Then, from eq. (1931) with the relation (19 7p . we find the deviation from the 4D SGUT 
prediction of the QCD coupling at M z , i.e. Aa s 



a KK _ a SGUT,0 as 



Aa s (M z ) 



~al(M z ){b\n-^ 

2l\ Mn_ 



B-L 



1 



In 



re" 2 




4 - 






U 


i 4 





1 ~ 6 

2 V 7 ; 



Q / l l. ,. 

7/i \iu)e 

iV 2 2 1 ' 



— nu/4 



4 i 



A 



(9£ 



The first term corresponds to the contribution due to the extra particles above the B — L 
scale. The second term is the volume dependent correction due to the KK massive modes 
while the third part containing the theta functions is the shape dependent correction. The 
last term A is the effect of the brane-localized gauge couplings. 

When u ~ 1, the shape dependent term is subdominant compared to the other log- 
arithmic terms. As can be shown explicitly in the specific models, the last term can be 
also ignored by making a strong coupling assumption at the cutoff scale. Then, the first 
two logarithms become a dominant contribution. For b(b — |) > 0, we can see that the 
individual logarithm can be large, being compatible with the gauge coupling unification 
due to a cancellation. 
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Now we consider the behavior of our result in the 5D limit with u = Rq/R^ 3> 1. In 
this case, the bulk gauge group becomes the Pati-Salam and there remain only two fixed 
points with the Pati-Salam group and the SM gauge group enlarged with a U(l) factor, 
respectively. Some relevant discussion on this limit has been made in Ref . [24] , concerning 
the power-like threshold corrections. Since |i?i(;z|Mi)| ~ 2e _7rU//4 | sin(7T2;)| for u ^> 1, the 
shape dependent terms could give a significant effect on the gauge coupling unification by 
the non-universal power-like threshold corrections. Thus, in this 5D limit, eq. fl98l) becomes 

Aa s (M z ) « -^al{M z )\l\n^--{b- Q -)HM,^V) 

-\(b - ln(4e^) + V+ M + A}. (99) 

Therefore, we can interpret that the effective 5D gauge coupling (l/pf = l/(glRe)) also 
gets a power-like threshold correction proportional to u/R 6 ~ I/R5 which is set by the 
mass scale of heavy gauge bosons belonging to SO(10)/ SU(4:) x SU(2) L x SU(2) R . 

On the other hand, when we take a different 5D limit for m Cl, the bulk gauge group 
becomes SU(h) x U(l)x and there remain only two fixed points with SU(5) x U (l)x and the 
SM gauge group enlarged with a U(l) factor, respectively. In this case, the appearing power 
threshold corrections are universal so there is no power threshold correction in eq. fl98l) . 

4.3 Gauge coupling unification in some 50(10) orbifold GUT 
models 

Now we are in a position to apply our general formula ( J98l) to particular cases for the 
unification of the SM gauge couplings. To this purpose, we consider some known SO (10) 
models of embedding the MSSM into the extra dimensions. In the minimal model(: model 
I) [13] that contains Higgs fields in the bulk for breaking U(l) b-l and the SM gauge groupie, 
there are 4 10's with parities (771,772) such as H\ = (+, +), H 2 = (+, — ), H 3 = (— , +) and 
if 4 = (—,—), and one pair of 16 and 16 with parities <3> = (—,+), $ c = (—,+). Then, 
the resulting massless modes are two doublet Higgs fields iff and if 2 from ifi and if 2, 
and Gl,Gi, (D c , N c ), (D, N) from H 3 ,H^,^ and $ c in order. Moreover, each family of 
quarks and leptons is introduced as a 16 being localized at the fixed point without SO (10) 
gauge symmetry. After the B — L breaking via the bulk 16's with (N) = (N c ) ^ 0, 
neutrino masses are generated at the fixed points by a usual see-saw mechanism. Moreover, 
G3, G4, (D c , N c ), (D, N) can acquire masses of order the B — L breaking scale by the brane 
superpotential [13, 14] W = XNDG C 3 + X'N C D C G 4 for (N) = (N c ) ^ 0. In this case, since 
Eio^i 10 = 0. Eie^i 16 = Ei6^i% 16 = -2, we get the values b = f , &"+ = b~ = f in 
eq. (|98p . Thus, in the 5D limit, because b~ + is nonzero in eq. fl99l) . there exists an effective 

6 In order to cancel the bulk anomalies due to one 45, we need to add in the bulk two 10's. So, it is 
necessary to have two Higgs doublets of the 10's in the bulk unlike in 5D case [3]. Moreover, in order to 
break the U(1)b-l, we need one 16 in the bulk. However, for cancellation of localized and bulk anomalies, 
one needs one 16 and two more 10's. 
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Figure 1: The la and 2a band of Aa s : the model I on the left and the model II on the 
right for u = Rq/R^ ~ 1. The dashed lines and the continuous lines denote la and 2a 
bounds of the experimental data, respectively. Only in the region above the straight line 
at M c /Mb-l = 1, the B — L breaking brane mass terms can be neglected. 



5D power-like threshold correction to the QCD coupling so the threshold correction is 
sensitive to the shape modulus. 

We consider another 6D 5*0(10) GUT model where the realistic flavor structure of the 
SM was discussed(: model II) [14]. In this case, on top of the minimal model, there are 
more hyper multiplets: 2 10's such as H 5 = (— ,+) and H 6 = (—,—), and one pair of 
16 and 16 with <fi = (+, +) and <p c = (+,+). Then, there are additional zero modes 
GI,Gq, L, L c from H^,HQ,(j) and C in order. They are assumed to get brane masses 
of order the GUT scale. Thus, the running of gauge couplings between the GUT scale 
and the B — L breaking scale is the same as in the minimal model. In this case, since 

Y^xoVl® = -2, Eie^i 16 = Eie^iV 16 = °> we S et the values ~ h = f > b ~ + = and b ~~ = f 
in eq. (158"j) . Thus, in the 5D limit, because b~ + = in eq. (1551) . there is no effective 5D 
power-like threshold correction to the QCD coupling. 

From the data of the electroweak gauge couplings at the scale of the Z mass, one 
can compare the predicted value of the QCD coupling in a theory to a measure one [25] 
a exp _ o.H76±0.0020. In the 4D supersymmetric GUTs, the prediction without threshold 
corrections for the QCD coupling is af GC/T, ° = 0.130 ± 0.004. Thus, in this case, there is a 
discrepancy from the experimental data as Sa s = cuf GUT >° _ a e s x P = 0.0124 ± 0.0045. 

First we consider the case with isotropic compactification of the extra dimensions, 
u ~ 1. In both models, since b = y, we can see from eq. (198|) that logarithmic contributions 
of zero modes and those of KK massive modes appear with opposite signs so that there is 
a possibility of having the large volume of extra dimensions consistent with perturbativity 
and gauge coupling unification. Ignoring the unknown brane-localized gauge couplings and 
the B — L breaking effect, we depict in Fig. 1 the parameter space of (M c , Mb-l) with 
u ~ 1, being compatible with the experimental data. If we take M*/M c ~ 63/a/C* ~ 
22 for strong coupling assumption^ at the 6D fundamental scale [8], the correction due 

7 We included the group theory factor C = 8 for the SO(10) bulk gauge group in the naive dimensional 
analysis compared to [8]. 
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Figure 2: The la and 2a band of Aa s : the model I on the left and the model II on 
the right for M*/M c ~ 22. The dashed lines and the continuous lines denote la and 2a 
bounds of the experimental data, respectively. Only in the region above the straight line 
at M c /Mb-l = 1, the B — L breaking brane mass terms can be neglected. 



to the brane-localized gauge couplings becomes A = 0(1) so it is negligible to the KK 
threshold corrections which is of order h^M*/ 'M c ) ~ 3. In the model I(II), for M*/M c ~ 22, 
Mb-l/M c can be as small as 0.23(0.12) at the 2a level. For M B _ L /M C < 1, the KK 
massive modes of the color triplets are modified to rri^ B n6 ~ (n 5 /2R 5 ) 2 + (n 6 /2R e ) 2 + 
cM B _ L where c is of order unity independent of the KK level for R 5 ^ i? 6 [26]. In 
this case, the B — L breaking effect to the differential running is estimated as M B _ L /M 2 
in comparison to A in eq. (198]) . so it is also suppressed compared to the KK threshold 
corrections. Apart from the two models, we can consider other possibilities of embedding 
the matter representations into extra dimensions, like in the field-theory limit of a successful 
string orbifold compactification [16] where there are two families at the fixed points and 
one family in the bulk. In view of the general formula fl98l) . however, as far as an extra 
particle contributes to the running of the gauge couplings above the B — L breaking scale, 
Mb-l tends to be close to M c for the success of the gauge coupling unification. 

Next we consider the shape dependence of the loop corrections. As in the previous case, 
we make the strong coupling assumption and take the B — L breaking scale to be below the 
compactification scale. Then, we depict in Fig. 2 the parameter space of (u, M B -l) with 
M^jM c ~ 22, being compatible with the experimental data. Therefore, we can see the 
clear difference between the two models, through the dependence of the compatible B — L 
breaking scale on the shape modulus. As shown in the limit of anisotropic compactification 
of the extra dimensions in eq. (1991) . in the model I, the power-like threshold correction in 
the effective 5D theory gives a sizable contribution to the differential running of the gauge 
couplings, thus the B — L breaking scale is more sensitive to the shape modulus. 

5 Conclusion 

We have considered the one-loop effective action for gauge bosons in six-dimensional orb- 
ifold GUT models with a number of hyper multiplets satisfying arbitrary local discrete 
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twists. From the obtained effective action, we encountered the divergences which require 
the introduction of brane-localized gauge kinetic terms and a bulk higher derivative term. 
Moreover, we derived the general expressions for the running of the gauge couplings of 
zero-mode gauge bosons in the low 4D momentum limit. Since the KK massive mode cor- 
rections depend on the parity actions, in general they can give a non-universal contribution 
to the gauge coupling running. 

By taking a concrete example such as the 5*0(10) orbifold GUTs, we estimated the 
corresponding KK massive mode corrections to the QCD coupling at the scale of the Z 
mass. The extra U(l) or the B — L symmetry is broken by the VEV of bulk singlets below 
the unification scale. Then, the extra color triplets, which accompany bulk singlets for 
a full 16, also appear as zero modes so that they can lead to an additional logarithmic 
running of the gauge couplings starting at the B — L breaking scale. 

In the case with isotropic compactification of the extra dimensions, there is a partial 
cancellation between two dominant logarithmic corrections; the volume dependent part of 
the KK threshold corrections and the threshold corrections due to the extra color triplets. 
In this case, we argued that the large volume of the extra dimensions can be compatible 
with the gauge coupling unification and perturbativity. We also considered the case with 
anisotropic compactification of the extra dimensions for which a 5D orbifold GUT limit 
can be discussed. In this case, the shape dependent part of the KK threshold corrections 
corresponds to non-universal power-like corrections in the compactification scales. These 
power-like corrections are calculable because they are finite in the 6D sense. This situation 
is in contrast to the genuine 5D orbifold GUTs with non-simple groups where power-like 
corrections with the cutoff dependence is uncalculable. Further, we showed that for a fixed 
volome of the extra dimensions compatible with a strong coupling assumption, the B — L 
breaking scale is sensitive to the change of the shape modulus, in order to maintain the 
success of the gauge coupling unification. 
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Appendix A: Propagators on GUT orbifolds 

Suppose that the bulk gauge bosons satisfy the boundary conditions, 



ViA^z)Vr = PiA„{-z + Zi)P: 
ViA m (z)Vr l = -P % A m {-z + z i )P i 



-i 



-i 



A^z + Zi), 
A m (z + Zi) 



(A.l) 
(A.2) 



where Pf — l(z = 0, 1, 2, 3) and [Pi, P 3 ] = 0. Then, we can write the gauge bosons Am on 
orbifolds in terms of gauge fields Am having AirR^^ periodicities along the extra dimensions 
as 



A M (z) 



n 



;(1 + V, 



-In 



(A.3) 



After taking into account the relation (fTTI) . the field redefinition becomes in terms of 
component fields in the group space 



S a h Al{z) + ^(3,)" b A%-z + 2zi) + ^ (QiQi) a + 2z % - 2^- 



i=0 



+ (QoQiQ2) a b A b J-z + 2z - 2 Zl + 2z 2 ) 



(A.4) 



5\A b (z) 



i=0 
\a /i 6 



£(Q0 8 *C(-* + 2*) + E (m-)%<(^ + 2^-2^) 



(Q QlQ 2 ) a + 2^0 - 22 X + 2Z 2 ) 



(A.5) 



where (Qi) a b = tr(T a PjT;,Pj). Thus, the orbifold-compatible functional differentiations for 
gauge fields are 



\ d 12 ) b 



5 a b 5\ Zl - z 2 ) ± E(Q,) Q 6 <5 2 (>i + z 2 - 2zi) 



i=0 



±(Q QiQ 2 ) a b 5 2 ( Zl + z 2 - 2z + 2z x - 2z 2 ) 



■8\x 1 -x 2 ). (A.6) 



Consequently, the propagator of gauge fields in the Feynman gauge is given by 



(A a M ( Zl )A bN (z 2 )) = gMN(5t 3 M ) a cG(z 3 -z i )(5i N T b 



(A.7) 
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where G(z 3 — z±) is a bulk scalar propagator satisfying 4:irR 5 ^ periodicities. Then, we 
obtain the propagator of gauge fields in 6D momentum space as 

(Al( Zl )A bv (z 2 )) -> g^{l + Q Q lC os{2p 5 7rR 5 )) a C {1 + Q Q 2 cos{2p 6 7iR 6 )) c d 

j (5fip> + QpSf-f) d b _ jf ff\\a (\ o\ 

9 ^2 _ =9iw{l*g,+ {P,P,P)) b, (A.8J 

Zi p p 



{A a m ( Zl )A bn (z 2 )) -> ^ ro „^(l + QoQiCos(2p 5 7ri2 5 )) c (l + g Q2Cos(2p 6 7ri2 6 )) c d 



x- 



(Jipj? QoSp-p 1 ) b _ (r / -,fxx a /a n\ 
2 ^ = gmn{Lrg,-{P,P,P )) b (A.yj 



where p = (p^pe) = (n 5 /(2R 5 ),n 6 /(2R G )) with n 5 ,n 6 integers. We note that the Wilson 
line effects encoded into Q\ and Q 2 are factorized as two matrices in front of the propagator 
on orbifolds without Wilson lines. 

Next let us consider a complex scalar field in the fundamental representation, satisfying 
orbifold boundary conditions, 

Vi<l>(z) = ruPi^-z + z^ = 4>(z + Zi) (A. 10) 

with r]i = +1 or —1. Similarly, we can write the complex scalar field on orbifolds in 
terms of a complex scalar field satisfying 47ri? 5i6 periodicities as 



4>{z). 



(A.11) 



Then, we can write the above equation in terms of component fields as 

1 r 2 

r(z) = - 6 a b <Kz)+Y,Vi(Pi) a b<l>(-z + 2z i )+ J2 ViV j (PiP j ) a b^(z + 2z i -2z j ) 
L ?;=o i<ri^a 



+mmm(PoPiP2) a b<t>\-z + 2z - 2 Zl + 2* 2 ) 



(A.12) 



Therefore, the orbifold-compatible functional differentiation for a complex scalar field is 

2 



; 12) a 6 



5 a b 5\ Zl -z 2 ) + Y, Vi(Pi) a bS 2 ( Zl + z 2 - 2zi) 



i=Q 



+ ViVAPPjT bS 2 (z l -z 2 - 2zi + 2 Zj ) 

+VoViV2(PoPiP2) a b S 2 (z 1 + z 2 - 2z + 2z x - 2z 2 ) 
Consequently, the propagator of a complex scalar is given by 

$ a (*i)? 6 te)> = (sf 3 ) a cG(z 3 - z 4 )(siy b 



6\x 1 -x 2 ). (A.13) 



(A.14) 
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or in 6D momentum space, 



^ a { Zl )^ b {z 2 )) -> ^(1 + r ]0 r ]l P P 1 cos(2p 5 7iR 5 )) a c (l + r ]0 r ]2 P P 2 cos(2p (i 7rR 6 )) c d 



i (s P ^+voPoS,,- f y b s ((V(Pi ^ ))v (A . 15) 

z p p 



For scalar fields of other representations, we only have to replace the parity matrices with 
the ones for corresponding representations. 

Finally let us consider a bulk left-handed fermion in the fundamental representation, 
satisfying the boundary conditions, 

VMz) = i£a 5 PM-z + Zi) = + Zi) (A.16) 

with £j = +1 or — 1. Following the similar procedure, the propagator of a bulk fermion is 
given in 6D momentum space as 

(r(z 1 )^ b (z 2 )) - i(l + ^i^oA cos(2p 5 7rP 5 )) a c (l + ^ 2 P P 2 cos(2p 6 7rP 6 )) c d 

if <5 d b 6 Pti ? d 8p- f 
X o 1 J~: ; ?o(-roJ ; ^75 



2U + 75P5+P6 ^ + 75^5+^6 

EE(^(p,p,^))° 6 . (A.17) 

On the other hand, for a bulk right-handed gaugino in the adjoint representation, the 
propagator in 6D momentum space takes the above form with £jPj replaced by Qi and 

^ + 75P5 + -> ^ + 75^5 - P6- 

Appendix B: KK summations in 6D orbifolds 

We consider the following KK summation (with c > 0, ai j2 > 0, < c± j2 < 1): 

-6/2 

7r[c + ai(ni + CiY + a 2 (n 2 + c 2 y\ 

ni,n2GZ 



Jo[c;ci,c 2 ] ee r[e/2] ^ [tt[c + ai(m + c x ) 2 + a 2 (n 2 + c 2 ) 2 ] 



/ -Kt\c+a 1 (n 1 +c 1 ) 2 +a 2 (n2+c 2 ) 2 } I R i \ 

/ t l- e/2 • ^- A J 

ni,n 2 ez u 



If < c/ai < 1, with notations 7(711) = z{n\)j yfa 2 — 7 c 2 ; and -2(711) = c+ai(77i+ci) 2 , 
u ee ^/ai/a 2 , s Sl ee 27rn 1A /c/ai, 7^ = 0.577216..., we obtain [6] (in the text ai = 1/P 2 , 
a 2 = l/P 2 ) 
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Jo[c] Ci,C 2 ] 
-$>|l-< 

ni6Z 



TTC 



-2 



■ + ln 



-27T7(ni) 



p>l 



An ai e Ts + ^( Cl ) + ^( _Cl ) 

I>+l/2] f-C 

(p+1)! L°i 



+ 27TM 



+ c\ - (c/ai + c?) : 



p+i 



C[2p+l,l+ Cl ]+C[2p+l,l- Cl ]) (B.2) 



while if we have c/ai > 1, then 
Jo[c; ci,c 2 ] 



TTC 



"-2 , 
hln 

e 




-J^ln l_ e -2-7(m) 






niSZ 



ni>0 



cos(27rni ci) 
hi 



Ms 



(B.3) 



The pole structure is the same for both cases; if c/a\ > 1 and except the first square 
bracket, no power-like terms in c are present (the last one being suppressed due to K\). 
Here Ct- 2 ? ] is the Hurwitz Zeta function and ip(x) = d/dx lnr[x] and K\ is the modified 
Bessel function. 

Finally, we quote here a limiting case for the behaviour of the function J Q 



Jo[c< 1;0,0] 



Jo[c< 1; 0,1/2] = 



TTC 



-2 



+ In 



Ane lE ai\q(i^ai/a 2 )\ 



-In 

7TC 



An 2 \r](i vV/g^)! 4 a 2 1 



— In c 



v /ala 2 " 



— + ln(47re^ B ai) - nu + 2 ^(-l) n In |1 - e 



—2irun\2 



n>l 



-In 



#i(l/2|iu) 



Jo[c< 1; 1/2,0] = 



7TC 



rj{iu) 



— + ln(7re _7B ai) - 2 In 



n>l 



1 + e- 



1-e" 



2n 



-In 



7?i(-iu/2|m) 



ITU /A 



rj{iu) 



Jo[c< 1; 1/2, 1/2] 



7TC 



\/ a l 2 

-In 



— + ln(7re- 7B ai) - 2 ^(-l) n In 



^(l/2-m/2|m) c _ TO/4 
rj{iu) 



n>l 
2 



1 


j g—nun 


2- 


1 


g — 7run 





(B.4) 
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Appendix C: Definition of special functions 

The Hurwitz Zeta function Ct- 2 ? ] is defined as 

CM = J2(n + a)- z 



n>0 



with Rez > 1 and a^O, —1, —2, 

The modified Bessel function K n is defined through 



I 



OO r\ 

v-l-bxP-ax-P Z 



dxx v t 



P 



2p 



with 



K x [x\ = e 



7T 




2x 



1 + — - — """"o + 0(l/x 3 ) 
8x 128x 2 V 1 ' 



(C.l) 



K*(2Vab), Re(b), Re(a) > (C.2) 



(C.3) 



which is strongly suppressed at large argument. 
In the text, we used the Dedekind Eta function 



2iiTT n 



n>l 

rji-l/r) = V^ V (r), V (r+l) = e^ 12 V (r), 
and the Jacobi Theta function i?i 

0i(*|t) = 2 ? 1 / 8 sin(7rz) -g n )(l -gV^)(l -g n e" 

n>l 

_ _■ ^2^_^n e iTT T (n+l/2) 2 e (2n+l)iTvz 



-2iirz\ 



n€Z 

which has the properties 



•& x {z\t + 1) = e i7r/4 tfi(z|r), 
+ ljr) = -0i(2|t), 



i?i(z + t|t) = 



-inr— 2inz 



A(~z/r\ - 1/r) = e^r^e^^Aizlr). 



(C.4) 



q = e 



2ilTT 



(C.5) 



(C.6) 



Appendix D: Some group theory for 50(10) GUT 
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• Relations between fundamental and other representations: 

For SU(N) and SO(2N) gauge groups considered in the paper, we have [20] 



trAdjFs£/(jv) - 2N tr N F SU ( N) , (D.l) 

tr a^Fsu(N) = ( N - 2 ) tr N Fg U( ^ N) , (D.2) 

tr^Fj^ = ±(N-2)(N-3)tY N Fl u{N) , (D.3) 

tr A dj-Fjo(2Ar) = 2(N - l)tr 2JV F| 0(2JV) , (D.4) 

tr 2 ]v-iFj o(2A — 2^ 4 tr 2 ArF| 0(2Ar) (D.5) 



where the subindex of the trace implies the representation of the group, for in- 
stance, a % ^(a l ^ k ) is the second(third) rank totally antisymmetric tensor representation 
of SU(N). In the text, we take the normalization, tr N (T a T b ) = \5 ab for SU(N) and 
tr 2N (T a T b ) = 5 ab for SO(2N). 

• Computation of traces: 

By standard representation theory, we can do the decomposition of the quadratic 
Casimir for an adjoint representation of SO (10): under the Pati-Salam, 

trAdj-Fjo(io) = tr A dji 71 | [ /( 4 ) + AtT 6 Fg V ^ + 12tr 2 F su{2) L + 12tr 2 F| t7(2)fl 

+ tr Adj F su(2) L + tr Adj F su(2) R (D.6) 

and under the Georgi-Glashow, 

trAdj-Fjo(lO) = tr Adj-^|t/(5) + tr W F SU(5) + tr W F SU(5) + % F U(1) X - ( D ' 7 ) 

Using the definition tr Ad j(QiFj o(10) ) = facdfbefV cd Qf F so(io) F so(w) and the fact that 
Qi is equal to +1(— 1) for Z 2 -even(odd) modes of gauge fields, we get 

trAdj(<2i-Ffo(io)) = tr Adj-^|{/(4) — ^ T e F su(4) ~ 12tr 2 F| (7(2)i - 12tr 2 F su{2) R 

+^Ad] F sU(2) L + ^Ad} F SU(2) R 
= -8tl 2 F SU{2) L - & tr 2 F SU(2) R ( D - 8 ) 

and 

^Adj(Q2 F SO(W)) = tr Adj F iu(5) ~ tr W F SU(5) ~ tr W F SU(5) ~ % F U(1) X 

= ^ F su (5 ) ~ 8^(i )x (D.9) 

and likewise 

trAdj(Q3-Fso(io)) = 4tr 5 F st7(5) , - 8F u{1)x/ . (D.10) 

Let us consider a similar decomposition of the index of the other representation of 
SO (10). First, the index of a fundamental representation of SO (10) is decomposed 
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as 



tr io-^lo(io) - tr 6 F su(4) + 2tr i F su(2) L + 2tr 2 F| (7(2)ij 

= fr5-F|tf(5) + 2^(l)x + tT 5 F SU(5) + 2 F U{l)x- ( D - n ) 

Then, with the parity matrices in the traces, we get 

tri (PiF| o(10) ) = -tTeFgy^ + 2tr 2 Ff {7(2)i + 2tr 2 F| (7(2)jj 

= -2tr 4 F| {7(4) + 2tr 2 F 2 77(2)i + 2tr 2 F| [/(2)j? (D.12) 

and 

tr 10 (P 2 F 5 , o{10) ) = -tr 5 F S{7(5) - ^ F u(i) x + tr 5 F su(5) + 2 F u(i) x = (D-13) 

and similarly tr 10 (P 3 Fg O , 1Q ,) = 0. Next, we also do the decomposition of the index 
of a 16 spinor representation of SO (10) as 

tr i6P|o(io) = 2tr 4 F| {/(4) + 2tr 4 F 5[/(4) + 4tr 2 F| [/(2)i + 4:tr 2 F su(2) R 

= trxo^B) + tr^^ + 1(10 + 45 + 25)F 2 (1)x . (D.14) 

Then, we get the necessary traces for a 16 as 

tri6(-PiP|o(io)) = 2tr 4 Fj [/(4) - 2tr 4 F| [/(4) + 4tr 2 F| [/(2)L - 4tr 2 F| [/ ( 2 ) fl 

= ^Fj^ - ^Fj^ (D.15) 



and 



tr 16 (P 2 F| o(10) ) = -tr 10 F| [/(5) +tr5Fj [/(5) + l(-10 + 45 + 25)F2 (lb 



and likewise 



= -2tr 5 i^ (5) + l F % {1)x (D.16) 

3 

^16(^3-^50(10)) = _2tr 5-P?l/(5)' + g^Wx'' ( D - 17 ) 

Relations between U(l) generators: 

There are three maximal subgroups of SO(10), Georgi-Glashow (377(5) x U(l)x) 
and Pati-Salam (517(4) x SU(2) L x SU(2) R ) and flipped S77(5) (SU(5)' x C/(l)x')- 
For a fundamental representation of 5*0(10), the U(l) generators are given by 



Y = diag(i, \,\,~\, ~\) x o 2 , * = diag(2, 2, 2, 2, 2) x a 2 , 

B — L — diag(| 1 1 0, 0) x a\ T 3R = diag(0, 0, 0, ~, - 1 -) x a 2 . (D.18) 
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Thus, we obtain the relation between £7(1) generators appearing in the different 
subgroups, 



Y = T 3R + ~(B-L), X = -4T 3R + 3(B-L). 



(D.19) 



Moreover, by comparing the flipped SU(5) to the Georgi-Glashow as N <-> e c L and 
u c L «-> d c L in 16 and h\ <-> /i2 in 10, we get another relation 



Y = -(-Y' + X'), X = i(24Y' + X'). 
5 5 



(D.20) 



Using the above relations, we can also find the relations between the U(l) gauge 
bosons. To this, let us consider the bulk kinetic terms for U(l) gauge bosons A Y , A x 
and a charged field 0: 



£bulk D 



L Fl - ±-F* + 



W Y w x 



IX 
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:d.2i: 



where g\ = gx at tree level. By writing 



\yA y + ^r=XA x 

5 



\y'A y . + -^X'A X , 
5 v/iO 



T 3R A R + J-(B-L)A B -l, 



(D.22) 



and using eqs. ( ID. 191) and ( ID. 201) . we obtain the relations between the U(l) gauge 
bosons as 



Ay 

A R 



~{-A Y + 2V6A x ), A X , = ^{2V6A Y + A X ), 



^{Ay-J^A x ), A b _ l 



A y + A X/ 



(D.23) 
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